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Abstract 

In this work an extension of the Gaussian model of the stochastic vacuum is presented. It 
consists of including higher cumulants than just the second one in the cluster expansion in QCD. 
The influence of nonabelian fourth cumulants on the potential of a static quark-antiquarkpair is 
examined and the formation of flux tubes between a static qq-p&ir is investigated. It is found that 
the fourth cumulants can contribute to chromomagnetic flux tubes. Furthermore the contribution 
of fourth cumulants to the total cross section of soft high energy hadron-hadron scattering is 
examined and it is found that the fourth cumulants do not change the general picture obtained in 
the Gaussian model. 



* This paper is dedicated to Franz Wegner on the occasion of his 60th birthday 



I. INTRODUCTION 



The linked cluster expansion |TJ , J2| is a powerful tool in statistical mechanics to evaluate 
expectation values. Especially if the measure for forming the average is not known it allows 
a systematic expansion. In particle physics this method is much less widely used. Here 
the weight, given by the exponential function of the action is known at least formally and 
therefore there are two principle ways to evaluate expectation values: 

The perturbative one where the action is split in a part quadratic in the fields, and the 
rest, which is expanded as a power series. The Gaussian integral due to the quadratic part 
can be performed and this yields the well known perturbation series 0,0]. 

Another way is to discretize the action in a gauge invariant way on the lattice || , || and 
perform the functional integration with numerical methods. 

In the model of the stochastic vacuum (MSV) [0, ||, |§ the linked cluster expansion of 
the long distance behaviour of QCD is the starting point. It turns out that if the stochastic 
variables are taken to be the gluon field strenghts a convergent cluster expansion yields 
automatically the salient feature of non-perturbative QCD namely confinement. This ap- 
proach has been successfully applied to many processes and it turned out that the simplest 
form of the cluster expansion namely the one where all cumulants higher than the second 
one are neglected yields a surprisingly good description of many phenomena ranging from 
charmonium spectroscopy to high energy scattering (see the reviews (H| , [ H ] , [|1^] , |13| and 
the literature quoted there). A process with vanishing higher cumulants is just a Gaussian 
one where the functional integrals can be performed explicitly and all can be reduced to one 
correlator. 

Though the Gaussian approximation is astonishingly successful, it is nevertheless worth- 
while to investigate which principally new effects can occur if higher cumulants are taken 
into account. This is the aim of the present paper. If one wants to study the gauge invariant 
field contribution of a static quark- ant iquark pair and high energy scattering of hadrons, 
one has to evaluate the vacuum expectation value of a product of traces of Wegner- Wilson 
loops. In that case the leading term is the expectation value of the product of four field 
strength tensors. A fourth cumulant can contribute directly to that expectation value and 
for that reason in this paper we concentrate on the fourth cumulant. 

It is structured in the following way: In the second section we shortly repeat some basic 
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features of the linked cluster expansion and some important technical features of the model 
of the stochastic vacuum (MSV). In section |T| we discuss some possible choices of higher 
cumulants and the way how we calculate the effects of these higher cumulants. We present 
the results obtained in this way and discuss them. 



II. LINKED CLUSTER EXPANSION AND THE MODEL OF THE STOCHASTIC 
VACUUM 

The quantity we are mainly interested in, is the Wegner- Wilson loop in QCD 0. 
It is defined as the vacuum expectation value of the gauge invariant path ordered closed 
line integral over the exponential of the color-potential 

(W[C}) = (Pe-^ib^W^) (1) 

where A jU (z) = J2A'j I \ a /2 is the SU(3)-Lie-Algebra valued gluon potential and P denotes 
path ordering. 

Using the non-Abelian Stokes theorem |14[], [[OJ, |IB| , we can express W[C] by a surface 
integral over the gluon field strength 

(W[C]> = (p se -*4/ s F ^M^M) ( 2 ) 

where the border of S is the closed loop C and da^ v is the usual oriented surface element of 
S at x; we choose in the sums always [i < v. The path ordering of equation ([!]) becomes here 
a surface ordering implied by the non-Abelian Stokes theorem. A possible surface ordering 
is given in fig. [I] . If we say in the following x\ < X2 ■ ■ ■ < x n we mean ordered according 
to surface ordering. In an Abelian theory (electrodynamics) the field strength tensor is a 
matrix of C-numbers and no path ordering is necessary. Then (^) is just the familiar Stokes 
theorem in 4 dimensions. 

The linked clusters are now defined by: 

<W[C]> = P 5 exp[-E2& 

)» (3) 

Lorentz and gauge invariance imply ((F All/ (x) )) = and the series starts with n = 2, the 
"Gaussian" cluster. 
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FIG. 1: Transition of path ordering to surface ordering with 
the non-abelian Stokes theorem. 

By expanding the exponentials and rearranging the terms one can show that the cumu- 
lants ((.)) can be expressed through the expectation values (.) 0. 

((F 1 F 2 )) = (F 1 F 2 ) ((F^Fs)) = (F 1 F 2 F 3 ) (4) 

and for xi < x 2 . . . < x n 

((F 1 F 2 F 3 F A )) = (F l F 2 F 3 F 4 ) - {F,F 2 ) (F 3 F 4 ) - (F 1 F 3 ) (F 2 F 4 ) - (F 1 F 4 ) (F 2 F 3 ) 

i (5) 

where F< = F^fa). 

For a Gaussian process all cumulants higher than the second one vanish by definition and 
we obtain directly for such a process 



(W[C]> 
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from which the area law of the W-loop can be derived easily |J, @ if we assume that the 
correlator itself falls off with a characteristic correlation length a. 

In order to evaluate expectation values of two (or more) Wegner- Wilson loops the simple 
factorization of the matrices does not suffice. In all applications studied so far, the following 
factorization hypothesis has been made for the color components of the field 



(F^F^F^F^) = (F^Ff) (F^F^) 

+ (F^Fg 3 ) (F 2 a2 F 4 a4 ) + (FfFf) (F 2 a2 F 3 a3 ) 



(7) 



with 

F? = F% Vi (x h w) (8) 
F^ iUi (xi,w) is the color-parallel transported field: 

F^x^w) = ^~ 1 (x i ,w)F IM/ (x i )^(x i ,w) (9) 

with 



w) = P exp 



ig / dX A M (w + X(xi — w)) (xf — w^) 



o 



(10) 



being the color-parallel transporter from x to w in the adjoint representation. Transforma- 
tion ([|) is necessary in order to give the non-local correlators a gauge invariant meaning. 

In the applications of the MSV to the calculations of more than one W-loop the following 
approximations are made: In the correlator (If^ vi (x^w)F^ U2 (x2,w)j) the dependence 
on the reference point w is neglected which is certainly only justifiable for a "reasonable" 
choice of the point w. A more technical assumption is the following: if the exponentials 
are expanded as power series in the correlator ((f^ i/i (xi,w)F^ 1/2 (x2,w) \\ only the term 
leading in (a/L) is taken into account, where L is the linear dimension of the loop and a the 
correlation length of the correlator. For a discussion and possible resolution of this rather 
ugly assumption see (JT7J, |T3 . 



Under these assumptions the following problems necessitating the evaluation of two loops 
have been solved: The calculation of the gauge invariant color field density J2 a ^(x) 2 , 
Y^ a B%(x) 2 [lTj and of total or diffractive cross sections in soft high energy reactions ||T8| . 
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In the next section we investigate in which way the very consistent results of the Gaussian 
model could be modified by the presence of higher cumulants. It is not astonishing that in 
that case the program cannot be performed as far as in the Gaussian approximation. We 
therefore concentrate on the leading effects of some higher cumulants. 

III. HIGHER CUMULANTS 
A. Choice of the cumulants 



In order to relax the factorization hypothesis of the MSV which only takes into account 
cumulants of order two, we introduce higher cumulants. The cumulant expression for a 
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correlator of fourth order is given by: 



fiai 77-02 pa 3 T?a A \ _ / / pai ™ 2 \ \ / / 77-03 77.04 

x /^l^l 1 ^2^2 A*3^3 W4 / \\ Ml^l M2^2 // \ \ W3 ^4^4 

— I — / / F a l P a 3 \\ // pa 2 77104 

1 \ \ W^l M3^3/ / \ \ ^2^2 ^4^4 



1 / / pai pa4 \\ // F a 2 pa 3 
1 W" 1 ^l^l M4^4// \\ ^2^2 M3^3 

+ //F ai F a2 F a3 F a4 \\. (11) 



We are not striving to find the most general ansatz for the fourth cumulant since a full 
inclusion of all fourth cumulants compatible with the general requirements of Lorentz and 
gauge covariance would introduce many new unknown parameters into the model. Even 



for the local condensate there are six independent ones pTJ • Any quantitative computation 
would lose its predictive power. Therefore we restrict ourselvers to studying the influence of 
some possible cumulants in a qualitative way. Of course the cumulants are not independent 



of each other, but there are some relations between them due to the equations of motion [21 



22| . However these relations are not sufficient to determine the higher cumulants through 
the lower ones. 

Two non-abelian example structures for a fourth cumulant, which satisfy the requirements 
of gauge and Lorentz invariance and which satisfy the normalisation condition: 

{{9**%^*%^)) = * {{9 2FF )) 2 . (12) 

are: 

' Vf;^ fa, w)f;i 2 ( X2 , w)f%j X3 , w)F^ 4 ( Xi , w) 



I 03 04 e { tfll V\ [13 U3 £(12 VI (14 ^4 6 Ml V\ (14 V4 ^(12 "2 M3 ^3 / 



and 
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^yf a\a 2 ef a 

~l _ /a3aie/a2a4e( e Mi^i/t4^4 e M2^2/i3^3 ^(11^1112^2^(13^3(14^4) 

~\~ f a 2 a z ef a 1 a i e{£(i 1 u 1 H2V 2 e (i3i J 3^4i'4 ~ ^l^iMS^A^i^/Mi^) } 

, Z 2 , Z3, Z4, Z5, Zq) (13) 



g 4 F^ Ui (x 1 ,w)F^ U2 (x 2 , w)F£ V3 (x 3 , w)F^ 4 ( Xi , w] 



■ V (( 9 >FF 
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+^0103^0204(^1^3^1^3 — ( ^A 1 1^3 ( ^M3^l)(^A t 2M4^2l y 4 _ ^ 7»2^4^ W^) 
+ ^104^0203(^1/44^1^4 ^/ii^4^4^i ) ^^2P.'i^V2Vi ^ fl2I^3^ ^3^2) } 

^(^1,^2,^3,^4,^5,^6)- (14) 

Here Zj are the differences between the coordinates Z\ — x\ — x 2 , Z2 = x\ — X3, . . . , z 6 = 
£3 — X4. The scalar function D falls off fast if any of the distances \zi\ becomes large, it 
is normalized D(0, 0, . . . , 0) = 1. We have extracted the dimensionful quantity (g 2 FF) 2 , 
the number which determines the strength of the full correlator ([11]). The first structure 
is similar to the four gluon vertex, f ai a 2 e being the structure constants of SU(3). The 
correlation lengths aj are defined by the following conditions: 

/ D{0,...,z i ,...,0)dz i = a i i = l... 6. (15) 
Jo 

In the following we put a« = a and all distances are taken to be in units of the correlation 
length a. In a Euclidean space we make a simple Gaussian ansatz for D: 

D( Zl , z 2 , z 3 , 34, z 5 , z 6 ) = e -i4+4+4+4+4+4)/^\ (16) 

where A = ^=a. As the distances z&, z^ and z e can be expressed as a function of Zi, z 2 and Z3, 
D is a function of the latter three independent variables only. Ansatz ([16| ) is not supposed 



to be a realistic choice but rather a simple manageable expression taking into account the 
cluster property. 

B. Single Wegner- Wilson loop: Potential 

In order to obtain the static potential of a quark- ant iquark pair, the expectation value of 
the Wegner- Wilson loop is calculated through the cumulant expansion 0,0]: 

(TrW[S]> = TrP s 

/ / da^^x^da^fa) ((F Mm (xi (x 2 ,w) )) 



cxp 



4-2! JsJs 

^ 2^4! Is Is Is Is da>MlUl ^ da ^ ( x ^) da ^ 3 (x3)da IHV4 (x 4 ) 

■ ( ( F Mm (xi , w)F fi2V2 (x 2 , w)F fi3U3 (x 3 , w)F fl4U4 (x 4 , w) ) ) 

9 6 
2 6 • 6! ' 



(17) 



where S is the integration surface enclosed by the loop and P5 denotes surface ordering. 

The first ansatz ([TJ) gives the following contribution to the potential, after performing 
the trace and summing over all color indices: 



V, 



(i) 



lim — Pc 

T^oo T 
1 



dcr^ lUl {xi)da t , 2U2 (x 2 )da lJL3U3 (x z )da^ UA (x i ) 



2 4 -4! Js Js JsJs 
1 

'T28 ' ' , -"'-/' : "' ! ' /'I ''J/'2"i!' 1 1 ■■'.>'■: 1 1 H-' ! 



JJ1 ILAM/i ^//.*i7yo//.QjyQ ^77.1 771 ll.nl/n^-tl.iIJi II.ATJA ~\~ 2 C7/1 7^1 //q 77^ 6 



J\f((g 2 FF^ D(zi,z 2 ,z 3 ,Z4,Z5,Ze). 



(18) 



Asuming that the Wegner- Wilson loop lies in a plane with the reference point also on this 
plane, the integration surface S 1 is just the rectangle enclosed by the loop. It is easy to see 
that the e-tensors vanish on this plane and that the cumulant of the first ansatz flllf ) does 
not contribute to the potential of the gg-pair. 

Inserting ansatz (|T4]) and summing over all color indices gives the following fourth cumu- 
lant contribution to the potential: 



V, 



(2) 



lim — Pc 



S JS JS JS 



dcr flll/1 (x 1 )da IM2l/2 (x 2 )da IMSI/;i (x 3 )da fl4U4 (x 4 ) 

" , "4 ~~ ^31/4^4^3) 



9216 ■ 2 4 ■ 4! 

1 16 fj, 2 ~ ^tliU2^(J. 2 iyi)(^fJ.3fl4^ 

~^{^fllfl3^Vlf3 ~~ ^Ml^^S^l) (^2M4^21'4 

+16(<5 /il/J4 5 UlU4 — 5^ jy 4 5^ Ul )( 5^ 2 ^ 3 5 U2 U3 — $ i_i 2 v?,$ HSV2) 

J\f ((/FF)) 2 D(zi, Z 2 , Zs, Z 4 , 26, *,)}• 



(19) 



After suitable parametrisation of the rectangular integration surface bounded by the loop, 



a non-vanishing contribution to the potential is obtained [£3 

5 



V, 



(2) 



lim 

t-^oo 3 • 2 25 iV, 



c 



JV (g 2 FF 



dn / dr 2 



dr-i 



dr A 



ds\ / <is 2 



ds* 



■T 3 .R 4 exp 



-R 2 



>2 - Si) 2 + (S3 - S 2 ) 2 + (S 4 - S3) 2 



—T 



+ (s 4 - s 2 ) 2 + (s 3 - s x ) 2 + (s 4 - sx) 2 
'r 2 - n) 2 + (r 3 - r 2 ) 2 + (r 4 - r 3 ) 2 



+ (r 4 - r 2 ) 2 + (r 3 - r x ) 2 + (r 4 - n) 2 



(20) 
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The integrations and the limit T — > oo can be performed numerically and the result is shown 
in figure |2|. The fourth cumulant contribution to the potential is plotted as a function of the 
quark- ant iquark distance R. For large R a linear rise is found, as in the Gaussian model. 
For short distances the potential behaves as i? 4 , in contrast to the Gaussian model, where 
the short distance behaviour is like R 2 . 




1 2 3 4 a 



FIG. 2: On the left the numerical result for the fourth cumulant contri- 
bution ( |l4|) to the potential of a static qq-p&ir is shown. For large 
quark separations R a linear rise is obtained whereas for small separa- 
tion a f? 4 behaviour is found. On the right the results of the Gaussian 
model for the Abelian (ab) and non-Abelian (n.ab) case are shown. 



C. The field of a static qq-pair 

We consider a static quark-antiquark pair in a color singlet state at fixed space points. 
It gives rise to a chromoelectric and chromomagnetic field whose components form the field 
strength tensor. In order to calculate the squared field strength we introduce a plaquette 
representing a small Wegner- Wilson loop in the /iz/-plane with center point x and 
dimension Rp (fig. |^). In the limit Rp — > its expansion gives the following relationship 
with the gauge invariant field density fll7| : 



P^Or) =N C - \RW £ W + °( R p) ( 21 ) 
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FIG. 3: Configuration of the Wegner- Wilson loop W and 
the plaquette P. The quark-antiquark axis is oriented in 
the a^-direction. The plaquette distance from the origin 
parallel to the quark-antiquark axis is represented by z and 
the distance perpendicular to the quark axis is x. In this 
figure the plaquette P34 is shown. 

where there is no summation over \x and v. This plaquette is evaluated in the presence of a 
static gg-pair whose world lines and color connectors form in turn a large rectangular loop 
of width R and length T going to infinity. We thus consider the quantity: 

4 (TrW^TrP^)) - (TrW[aS]> (TrP^(x)) 



fi 



X 



(TrW[dS\) 



(22) 
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(23) 



E and B being the chromoelectric and chromomagnetic field. The coordinates of the pla- 
quette are x, y and z: because of rotational symmetry in the Xi-avplane, we choose y = 0. 
Assuming that the Wegner- Wilson loop lies in the X3-X4-plane, x is the perpendicular dis- 
tance to the quark axis, whereas z represents the distance from the origin parallel to the 
quark axis see fig. [3[ We proceed by studying / M „ in the limit Rp — ► 0. After the expansion 
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of the exponentials in fl22|) we obtain: 
4 1 



/(*) 



Rpg 2 (TrW) 
, . ordered 



free 



E 



J 2"n! 



v — ^Tr 
2 2 2! 

ordered 



-E 



71=1 Z 



da - 



fiivi 



2 2 2! 



■Tr 



Tr[T ai ...T a 

rt" /T" 1 ! . . . F a ™ n 2 p a p b \ — 

i) M ii\V\ p n v n i) pu pa j 

.da™ v Tr[T ai . 
n n /? a i . . . F a " \ I n z F a F 

y - 1 /ill/1 - 1 Mn^u I \i> uv i 



d(T p,u d(7 pa 



fj,v pa i 

free 



d(T pu d(T pa 



b 

fJ,U^ pa 



(24) 



where the integrations over the surface bounded by the Wegner-Wilson loop are surface 
ordered. The matrix elements of / are given by the different orientations of the plaquette. 
In the Gaussian model it was possible to insert the factorizations and rewrite the summations 



in such a way as to perform the summations [0. Here this full summation is not feasible. 
Nevertheless to get an idea about what kind of contribution a fourth cumulant could give 
to the fields, we confine ourselves to the lowest non-trivial term. It should give a qualitative 



indication of the effects of higher cumulants. After inserting the cumulant definition (11) 
we obtain the leading contribution from the fourth cumulant: 

4 



5f :- 



R 4 P 9 2 
ordered 



free 



2 2 2! 
1 

222! 



pa( X p) 



Tr 



9 2 F^ 1 (x w )F^ U2 (x' w )g 2 F^(x P )F ( 



pa\ X P, 



(25) 



Since we are mainly interested in the x, z-dependence of Sf we omit the x and ^-independent 
normalization (TrW). 



To compute the contribution of each ansatz (|13|) or fll4|) the traces must be performed 
and a summation over all color indices must be carried out. After suitable parametrization 
of the plane surfaces bounded by the Wegner-Wilson loop and the plaquette, we obtain for 
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the contribution of the first ansatz 



31: 



5fi 



18432# 2 



and for the second ansatz (111) |23 



5f2 



-1 



147456# 



-M ((g 2 FF 



^0 10 0^ 

10 



yO 



( q 4 4 4\ 

4 4 4 

4 4 5 

4 4 5 



(26) 



(27) 



where both contributions have the same x, ^-dependence: 

/l rco 
R 2 dswds' w I d,T W dT' v 
1 J-oo 

■ exp 



'W 



Ax 2 + 2{z - -Rs w y + 2(z - -Rs' w ) 



+-R 2 (s' w - s w ) 2 + 2(t& + t#) + {j' w - r w f 



(28) 



The first ansatz ( |i~3"D only contributes to the chromomagnetic field in the direction parallel 
to the quark axis: 



5Bt 



N, 



c 



18432^ 2 



M((g 2 FF)) £{x,z,R) 



(29) 



whereas the second ansatz (14) contributes to all components of both the chromoelectric 
and the chromomagnetic field: 



SB 2 = SB 2 = 5Bt 



5E 2 = 5E 2 



1 -M((g 2 FF)) 2 i{x )Zl R) 



5E 2 



147456^ 2 
5 



N((g 2 FF))\{x )Z) R). 



(30) 



147456^2 

The numerical evaluation of the x, z-dependence £(x, z, R) of the fourth cumulant contri- 
butions to the fields is shown in figure ^ for the gg-separations R = 1 and R = 8. The result 
is that, for increasing separation between the gg-pair, a flux tube is formed. With (|29|) and 
(j30|) the following conclusions can be drawn: 
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FIG. 4: The function £(cc, z, R) in arbitrary units for the qq- 
separations R = 1 and R = 8. The quarks have coordinates 
x = and z = ±0.5a for i? = 1 and x = and z = ±4a for 
R = 8. One observes the formation of a flux tube, note the 
different z-axis scales of the two plots. 

• The first ansatz ( |T3|) leads to the formation of a purely chromomagnetic flux tube for 
the ^-component of the field, but to no area law of the Wegner- Wilson loop. 

• The contribution (^) of the second ansatz (|i~4[) can be split into a deformation of 
the vacuum field density with no preferred field components and a flux tube for the 
color-electric field in the z-direction. Of course both contribute to the energy density 
of the system. In this case confinement also reflects itself in the area law. 

We thus find for the first ansatz (|i~3|) a contribution to a chromomagnetic flux tube which 
was not present in the Gaussian approximation. The easiest interpretation of this behaviour 
is that the flux tube has been created by a pair of magnetic monopoles. At first sight 
there seems to be a contradiction between the formation of a flux tube and the absence of 
a contribution to the area law of the Wegner- Wilson loop. However in the Wegner- Wilson 
loop the exponent gA^dx 11 only takes into account the interaction of a color- electric charge. 
Therefore it is not astonishing that the area law of the Wegner- Wilson loop ([I]) is not a 
criterium for confinement of color- magnetic monopoles. 
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D. High Energy Scattering 



Soft hadron-hadron high energy scattering has been extensively studied using the Gaus- 
sian model fTPf , |TT| , flT2| , [ p"3|1 , POfl . We have investigated the effect of relaxing the Gaussian 



approximation in much the same way as for the static gg-pair. Both fourth cumulants ([13]) 



and (H) can be continued to Minkowski space in exactly the same way as the second cu- 
mulant and in the final results only the separations in transverse space enter; this means 
that the Euclidean expressions for the correlators can be used. In fig. [5] we give the total 
dipole-dipole cross section as a function of the dipole size for two equal dipoles. The quali- 
tative behaviour is the same as for the Gaussian model: at short distances the cross section 
is proportional to the fourth power of the size of the dipoles, for large sizes it approaches a 
quadratic dependence. We have also studied the dependence of the cross section on dipole 
size with one fixed-size dipole. Again we find the same behaviour as in the Gaussian model 
namely for small sizes a dependence on the square of dipole size and for large sizes an 
approach to a linear behaviour. 

IV. SUMMARY AND OUTLOOK 

We have found that the contributions of higher cumulants in the cluster expansion do not 
lead to a change of the qualitative behaviour of phenomenological results of the Gaussian 
model. There are however interesting new features in the field configurations of a static qq- 
pair. The resulting flux tube of ansatz ([14]) consists of an increased isotropic vacuum density 
inside the region described by the function £(x, z, R) and in addition of a color-electric flux 



tube. The ansatz (13) leads to the formation of a color- magnetic flux tube, very similar 
to the color-electric flux tube of the Gaussian model. The easiest interpretation of such a 
flux tube is that it has been generated by a pair of confined monopoles. Though a static 
bare quark carries only color-electric charge it is not impossible that it acquires a magnetic 
charge in the vacuum. For a color-magnetic monopole without color-electric charge the area 
law of the Wegner- Wilson loop as evaluated in the MSV is then no longer an indication of 
confinement. It would therefore be very interesting to investigate higher cumulants on the 
lattice since they might reveal new features of non-Abelian gauge theories. 
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[< g 2 FF > 2 a 10 Af] 



[< g 2 FF > 2 a 10 Af] 



1 2 3 4 5 a 



[< g 2 FF > 2 a w ] 




FIG. 5: The first two figures give the fourth cumulant contributions to 
the total dipole-dipole scattering cross section for ansatz (13) and ( |l~4D 
respectively. The last figure gives the results of the Gaussian model for 
the non-Abelian (n.ab) and the Abelian (ab) case. 
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